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We present a novel approximation for the treatment of two interacting magnetic impurities im- 
mersed into a noninteracting metallic host. The scheme is based on direct perturbation theory with 
respect to the hybridization between the impurity and band electrons. This two-impurity enhanced 
noncrossing approximation can fully incorporate the indirect interactions between the impurities 
which are mediated by the conduction electrons as well as additional arbitrary direct interaction 
matrix elements. We qualify the approximation by investigating the uncoupled case and conclude 
the two-impurity approximation to be equally accurate as its single impurity counterpart. The 
physical properties of the two-impurity Anderson model (TIAM) are first investigated in some lim- 
iting cases. In a regime, where each of the uncoupled two impurities would exhibit a pronounced 
Kondo effect, we ignore the indirect coupling via the conduction band and only incorporate direct 
interactions. For a ferromagnetic direct exchange coupling, the system displays a behavior similar 
to a spin-one Kondo effect, while an antiferromagnetic coupling competes with the Kondo effect and 
produces a pseudogap in the many-body Kondo resonance of the single-particle spectral function. 
Interestingly, a direct one-particle hopping also produces a pseudogap but additionally pronounced 
side-peaks emerge. This gap is characteristically different from the case with antiferromagnetic cou- 
pling, since it emerges as a consequence of distinct Kondo effects for the bonding and anti-bonding 
orbital, i.e. it reflects a splitting of even and odd parity states. For the general case of only indirect 
coupling via the conduction band the results show signatures of all the previously discussed limiting 
cases as function of the impurity-impurity distance. Oscillatory behavior in physical quantities is 
to be expected due to the generated Ruderman-Kittel-Kasuya-Yosida (RKKY)-interaction. We are 
led to the conclusion that the well known Doniach-scenario captures essential aspects this model, 
but the details especially at small distances are more complicated. 



PACS numbers: 71.27.-|-a,75.20.Hr,71.45.Gm,73.21.La 

I. INTRODUCTION 

The investigation of strong correlation effects plays a 
major role in contemporary condensed matter research^ 
as, for example, in the field of high-temperature cuprate 
superconductors,^ heavy-fermion systems'^ or frustrated 
magnetism.'* Due to the complexity of the phenomena 
and the theoretical models there arc still fundamental 
questions to be addressed, despite long and intensive re- 
search in these fields. 

Quite generally, the reason for this can be found in the 
presence of competing physical tendencies, whose char- 
acteristic energy scales are close together.^ Then, a very 
rich variety of physical scenarios results since subtle dif- 
ferences in the external conditions or small parameter 
changes can tip the balance from one dominating mech- 
anism to another. 

We consider the two-impurity Anderson model^ 
(TIAM) as one of the most basic models incorporating 
such competing interactions. There, two magnetic im- 
purities are immersed with finite distance into a host 
with a noninteracting conduction band and hybridiza- 
tions between conduction electrons and local interacting 
electrons. 

A single impurity is prone to the Kondo effect,^ where 



an effective antiferromagnetic exchange between the con- 
duction electrons and the local spin of the magnetic im- 
purity leads to a dynamic spin-screening. This archetyp- 
ical many-body effect is associated with an characteristic 
energy scale Tk (Kondo scale) which is nonanalytic and 
exponentially small in terms of the effective exchange. 

Additionally, the noninteracting band electrons in- 
duce an effective magnetic exchange between the two 
impurities,^ known as the Ruderman-Kittel-Kasuya- 
Yosida (RKKY)-interaction. This RKKY-exchange can 
be ferromagnetic or antiferromagnetic, depending on the 
spatial separation of the two impurities. 

The RKKY-exchange is expected to be dominant for 
a small hybridization between the impurity and the con- 
duction band and favors an inter-impurity coupling of 
the spins usually associated with long-ranged magnetic 
order. The Kondo coupling, on the other hand, strongly 
increases toward larger hybridization strengths and pro- 
motes individual Kondo effects at each impurity resulting 
in two uncoupled and screened impurities. This scenario 
proposed by Doniach^ for the competition between these 
two mechanisms is at the heart of the rich variety ob- 
served in many physical systems.'^ 

The TIAM as the simplest model for this competi- 
tion has be studied intensively-'^^"*'^ in the past and has 



2 



been realized in a recent experiment.^ The thermody- 
namic properties of the model are governed by a quan- 
tum critical point separating the ordered and paramag- 
netic (Kondo) ground states in case of identical impuri- 
ties and particle- hole symmetry. However, the quantum 
critical point is unstable and turns into a crossover once 
one of these symmetries is broken or a direct hopping is 
included. ^^^^'^ The dynamic properties under general con- 
ditions, that is varying impurity-impurity distance, hy- 
bridization and interaction strength, are not as clear. Es- 
pecially the situation with finite non-diagonal hybridiza- 
tion is not well studied since usually the RKKY exchange 
is simulated via a direct exchange interaction. 

In this work we present an approach to the a general 
TIAM which is formulated in the framework of direct 
perturbation theory in the hybridization^^ which utilizes 
skeleton graphs for families of time ordered diagrams. 
These techniques have been laid out in the literature, 
e.g. in connection with the Kondo effect or the mixed va- 
lence problem, and have served to define approximation 
schemes of great usefulness and quality.^^'^® 

There have already been approaches to the TIAM with 
these techniques, but these were restricted to infinite 
Coulomb interaction, [/ — >■ oo, neglected vertex correc- 
tions, and ignored non-diagonal ionic propagators (see 
below for details). The present scheme is an exten- 
sion of the finite- J7 enhanced noncrossing approxima- 
tion (ENCA)^°, which improves on the well known non- 
crossing approximation (NCA)^-'^ by the inclusion of low- 
order vertex corrections. It is also a direct extension of 
the ENCA for the multi-orbital situation as described in 
Rcf. 22. 

Apart from an application to the rich and interesting 
physics encountered in the TIAM, our new two-impurity 
solver also is of interest for extensions of the dynamical 
mean-field theory to include nonlocal two-site correla- 
tions. Compared to other existing two-impurity solvers 
the present scheme has some advantages as well as short- 
comings: Like all the schemes based on perturbation the- 
ory in the hybridization, it has some problems to cor- 
rectly describe the Fermi liquid properties at very low 
temperatures. This does not represent a major draw- 
back as many scenarios and questions of current inter- 
est can be addressed at temperatures accessible within 
these schemes (see, for example, Ref. 23). As will be 
demonstrated below and as it is known from the single- 
impurity casc,^^'^"'^^ these schemes correctly describe the 
non-perturbativc Kondo-physics and arc able to repro- 
duce the exponentially small low-energy Kondo scale. 
Additionally, there are no adjustable parameters and this 
theory is directly formulated for correlation functions at 
the real frequency axis This avoids the inaccuracies con- 
nected with a numerical analytic continuation of imagi- 
nary time data, as obtained, for example, from quantum 
Monte Carlo schemes. 



II. MODEL 

We consider two impurities (labeled as /-electrons) 
placed at at lattice sites Rj with j = {1,2}. That part of 
the Hamiltonian preserving all local occupation numbers 
is 

+ \ U[,nj^nl,. (1) 

The operators fj^ and /J^ are the usual annihilation and 
creation operators for impurity /-electrons with spin a = 
{+, -} = {t,;} at the lattice J and n^^ = /j^/,-^. The 
density-density Coulomb interaction is incorporated via 
the matrix elements Uj for on-site and for inter-site 
interaction. We also allow for more general couplings 
between the two impurities 

w^^Yl t^ifiJ^'^ - ^ E -^y^' 4 + ^'(^' (2) 

where the first term is direct single-particle hopping with 

amplitude tij and the second term a Heisenberg exchange 

-f 

with exchange coupling Jij and S_j the vector /-electron 
spin operator. The third term allows, e.g., correlated 
hopping (associated with matrix element X) and pair- 
hopping (associated with matrix element Y). 

These impurities are immersed into a lattice with a 
noninteracting conduction band (c-electrons) with the 
dispersion relation ej:., 

with the number operator for c-clcctrons n^^ = c|,^Cj,^ 
with spin a and crystal momentum k. Both subsystems 
hybridize via a term 

The total Hamiltonian then reads 

H = Hi + Wf + + V. (5) 

Even though in principal there is no limitation in ap- 
plying our calculational scheme to this model in full gen- 
erality, in this work we will mostly consider identical im- 
purities and focus on the half-filled particle-hole sym- 
metric case, where two electrons are placed in the two- 
impurity cluster. 

III. NOVEL TWO-IMPURITY SOLVER BASED 
ON DIRECT PERTURBATION THEORY 

The general idea of direct perturbation theory is to 
treat the two- impurity subsystem + H^^ + exactly 
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and to use V as the perturbation. With = V = 
the Fock-space of /-electrons is diagonahzed in terms of 
the local occupation numbers n^^. An eigenbasis of Hq 
is generated by all product states with the factors taken 
from the two versions of the local /-eigenbasis for the two 
sites j = 1 and 2 , respectively 
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The inclusion of the term makes it necessary to trans- 
form to a different many-electron eigenbasis of /-states 
in order to diagonalize the two-impurity subsystem. The 
following discussion of the novel two-impurity solver is 
formulated in general terms. It uses the afore mentioned 
multi-electron eigenstatcs of the diagonahzed /-system, 



hI\M)^Em\M), 



(7) 



labeled by quantum numbers M. The corresponding 
Hubbard-transfer operators 



Xmm- = \M){M'\ 



(8) 



allow for the decomposition of one-particle creation- and 
annihilation operators (with coefficients specified later), 
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(9) 



M.M' 



M,M' 
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In direct perturbation theory the partition function 
and the fermionic two-point Green functions are ex- 
pressed as contour integrals in the complex energy 
plane: 



dz 



M 



Tr 



7^(z) 



(10) 



-ImPji/;M(w + iS) 



1 f dz 

Ga.b(«w„) = ^ TT^e-'^^Tr n[z) A TZ{z + lUJn) B 

Z / ZTTl 



dz 



^ M',M", ^ 
M'" ,M 

{z, Z + iuJn)PM":M"' {z 4- ii^n)- 



(11) 



z-H 



and 



with the resolvent operator TZ{z) 

Zc = Trc ^e"^^''^ is the partition function of the non- 
interacting subsystem of band electrons. The Pm,M'{z) 



are so-called ionic propagators which capture the corre- 
lated dynamics of a transition from an ionic /-state with 
quantum number M' to a state with M . These are exact 
equations if the exact vertex functions Ajj/.jvj" (z, z + iujn) 
are used. In practice, however, the exact vertex func- 
tions are unknown and approximations for KmI:m" are 
employed. 

For simplicity, wc will confine our considerations to 
one-particle /-Green functions where A — fj^ and B = 

fj^,. For the two- impurity model the one-particle Green 
function is conveniently written as a 2 x 2 matrix for each 
spin component, 
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A central quantity of physical interest is the one-particle 
spectral function which is given by the diagonal part of 
the Green hmction. 
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ImGj? ^"t (i^ + i6), 



(13) 



where the limit < (5 — > is implied. 

In contrast to what is commonly needed for (effective) 
single-impurity models, ^""'^^ we introduced non-diagonal 
ionic propagators. 



P, 
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-{-f,M\n{z + E^)\M\-/) 



(14) 

which together with the diagonal components and the 
vertex functions AM;M'iz, z') are the central quantities 
of interest which need to be computed. 

The perturbational processes contributing to the ionic 
propagators and vertex functions take place along the 
imaginary time axis < t < /3 = ^. These processes 
can be partially summed-up to families of time-rotated 
diagrams on cylinders where the endpoints t = and 
T = /3 are identified. The Green function of Eq. (11) is 
represented by one particular diagram where the exter- 
nal line (accounting for the operator B) is attached to a 
bare vertex occuring at the latest time; the dressed ver- 
tex Ajv/ j\/'(z, z'), where the external line associated with 
the operator A is attached, occurs at an earlier time. The 
initial and final ionic state M is displayed twice in graph- 
ical representations for clarity, but is included only once 
in the corresponding formulas. 

As usual, propagators are favorably built up from irre- 
ducible self-energies I]jv/;a/' (z), which together with the 
vertex functions form a closed set of skeleton equations. 
Diagrammatic contributions to the self-energies are iden- 
tified as pieces of diagrams which cannot be separated 
into disjoint parts by cutting only the time axis of each 
site in the cluster. 

It is important to realize here, that the time-ordered 
processes, generated by expanding the resolvents in Eq. 
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(14) in terms of the hybridization V, generally extend 
over all sites of the cluster. Local pieces on one particu- 
lar site j can occur in proper time-order, either after an 
electron has is transferred from a different site I or after 
an interaction process, both induced by . If would 
have been regarded as a perturbation, too, these events 
would appear in the diagrams explicitly. However, we will 
take W-^ into account exactly via a non-diagonal form of 
the ionic propagators p''^^.j^p{z), which corresponds to 

the unperturbed Hamiltonian + Hq + . 

As a consequence of the all-embracing unique time- 
order in the cluster only one energy variable is present 
after the Laplace-transform. Another consequence of 
this fact concerns the dynamical interdependence of local 
parts of a process in a cluster. As already mentioned, all 
events occur in a definite time-order. Even if two pro- 
cesses happen to be located at separate impurities, e.g. 
the excitation of two band electrons and re-absorption at 
the respective same site, different time-orderings of the 
individual two excitations and absorptions are counted 
separately in this approach. For example, interchanging 
the emission times of the two processes generates or re- 
moves one crossing of two band-electron lines. As the 
approximation is organized in numbers of crossing band 
electron lines higher-order vertex corrections with a mul- 
titude of crossing lines are not included. This introduces 
artificial inter-site correlations. 

Therefore, all such approximations within a two- 
impurity model will in principle not faithfully reduce to 
two independent sites in the case of uncoupled impurities. 
They have to be thoroughly checked in the limiting case 
of two independent sites, where an ordinary Kondo effect 
is expected at each site. As it turns out, the two-impurity 
solvers to be introduced below are capable of describing 
this limiting case reasonably well. On the other hand, 
these approximation schemes will perform even better 
when real physical inter-site correlations make it neces- 
sary to respect overall time order. 

The connection between the irreducible self-energies 
and the ionic propagators is provided by a Dyson- 
equation 

Pa/;M'(^)=p1°Jm'(^) (15) 
+ Pm]m"^m";M"'{z)Pm"';M'{z)- 

AI",M"' 

As already stated above, the action of W-^^ is included in 
P^M-M' 1 i-C- it is a solution to the equation 

P^m]m-{^)^^t^5m,m' (16) 

Z - iiM 

A/,A/"^l/'';J\/'(-^)- 

^ ^''^ M" 

The irreducible self-energy Yim^m' incorporates the addi- 
tional modifications brought about by the hybridization 
V . The limit of an isolated cluster not hybridizing with 



the conduction band is thereby exactly reproduced for 
vanishing hybridization with the band states, since 

Vk~^Q: Ym-m'{z)^Q. (17) 

For the present case of a cluster with two sites j = 1,2 
and s-shcUs only, the quantum numbers are specified as 
M = (7711,7712) with iTij = {0,t,i, 2}. The annihilation 
operators of Eq. (9) are then given as 

/i.=[|0)ii(^|+a|-'T)ii(2|]®i2, (18) 

/2.=ii® [10)22(^1+^1-^)22(21], 

h =|0)„-(0| + (-1)^' \^)n{^\ + I2),,(2|. 

a 

The differing sign in the definition of the local identity 
operators for impurity 1 and 2 stems from the definition 
of the basis states, in particular from the prescribed order 
in which the creation operators act to produce the other 
basis state from the two-site vacuum. 

The overall sign of a diagram turns out to be more sub- 
tle for the two-impurity case than for a single-impurity. 
In the latter the overall sign is determined by simply in- 
cluding a factor of (—1) for each crossing of two band elec- 
tron lines. For the two-impurity case, we have found no 
better way than to count all exchange processes of Fermi- 
operators in the course of reducing expectation values of 
operator products to normal order. This is done algorith- 
mically on a computer along with the enumeration and 
classification of all different diagrams contribution to the 
quantity of interest. 

The analytic contributions are visualized in a diagram- 
matic language similar to the one developed long ago 
for the single-impurity problem, however, with a few 
characteristic differences. Even though there exists only 
one energy variable z for both sites, we distinguish two 
vertical lines corresponding to local processes on sites 1 
(left) and 2 (right) for clarity. These become decorated 
with vertices ordered along the vertical direction due to 
hybridization events. Bare vertices are depicted as cir- 
cles with the corresponding number of the site where the 
external line is attached at, and dressed vertices are de- 
picted as triangles; they embrace both lines as shown in 
Fig. 1(a) and (b). 

The ionic propagators between vertices are nonlocal 
objects and non-diagonal in quantum numbers M, and 
as such they describe the combined time evolution of both 
ionic states. Band-electron lines are always considered to 
be external to the whole cluster and can be tied to either 
of both sites indicated by the number inside the vertex 
symbol. The short horizontal lines denote amputation 
as usual, where the initial and final states are not in- 
cluded in the analytical representations. The translation 
of a diagram into an analytical expression is handled in 
complete analogy to the diagrammatic rules stated in the 
literature, ^^'^^ except for the sign rule mentioned above. 

Hybridization events occur as diagonal (j = /) and 
non-diagonal (j 7^ /) processes characterized by the (com- 
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(b) 
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FIG. 1. (a) Bare vertex depicting an absorption of a band 
electron with spin a at site 1, changing the two- impurity quan- 
tum number from M = (0, 0) to the M' = (a, 0). (b) Dressed 
vertex for the absorption of a band electron with spin a at 
site i. The left (right) vertical line denotes site 1 (site 2) and 
the incoming arrow a band electron. A dashed vertical line 
represents an empty and a full vertical line a singly occupied 
shell, where the spin of the electron is indicated next to the 
line. The number in the vertex symbol denotes the site, to 
which the band electron line is attached. For the bare vertex 
this can occur only on site 1, while for the dressed vertex this 
is not specified. In an analytic expression, the bare vertex 
amounts to a factor of 1, while the dressed vertex stands for 
the function A.qq.q„{z,z -\- e), when the band electron has an 
energy e and the local two-site cluster has the complex energy 
variable z. 



plcx) functions 



fe 



(19) 



Their imaginary part at the real frequency axis (z = 
w ± iS) 

A,- , iuj) = i [Fj- i - i5) - r,- , [uj + iS)] (20) 



is called hybridization function and is the physically most 
relevant part. The diagonal component at zero frequency 
corresponds to the well-known Anderson width of a res- 
onant impurity level, A = Ajj(a; = 0). As with the 
other quantities, the diagonal and non-diagonal compo- 
nents are conveniently treated as matrices. 

We refrain from presenting the complete system of self- 
energy and vertex-equations, which even after using all 
symmetries of the two-site Anderson-model, is far too 
large to present (It can be found, however, in Ref. 25.) 
Instead we exemplify the structure by selected examples 
for self-energy and vertex equations. 

The cluster vacuum with no electrons on each of the 
two ionic shells experiences excitations according to the 
self-energy diagrams shown in Fig. 2. The corresponding 



analytic contribution is 



(21) 



A,,2(e)/(e) A^o.o,(2, z + e) Paa:Oa{z + e) 
+ A,,i{e)f{e) Aj)o.o^(z,z + e) Poa-aoi'' + ^) 
+ A,,2(e)/(e) A'ao-fiaiz, z + e) PaQ-Mz + e) 
+ A,,i(e)/(e) K'{z, z + e) P^oMz + e) 



Among the graphical elements in Fig. 2 is an insertion 
with a diagonal line inside the square. This represents a 
non-diagonal ionic propagator where an electron is trans- 
ferred from one site of the cluster to the other. Fig- 
ure 3 presents the diagrammatic contributions to a non- 
diagonal self-energy. 

Up to now, all the equations for propagators and self- 
energies are exact. However, the dressed vertex func- 
tions, as indicated by the triangles in the graphs, have 
not been specified yet. Clearly, the exact vertex func- 
tions are in general unknown and approximations have 
to be introduced. 

The simplest approximation amounts to ignoring all 
vertex corrections and just replacing every full vertex by 



a bare one, that is Ai 



M-M' 



1. The result is the natural 



generalization of the NCA for finite U (usually termed 
SNCA) to the two-impurity case and will be referred to 
as two-impurity SNCA. 

In order to incorporate the correct order of magnitude 
for the Schricffer- Wolff exchange coupling for finite C/, the 
first order vertex corrections with one crossing of band 
electron lines have to be included. The extension of 
this ENCA for the single-impurity to the present two- 
impurity set-up will be termed two-impurity ENCA. As 
an example we present the vertex function Kqq.^q{z, z') 
within the two-impurity ENCA in Fig. 4. 

A closed system of coupled equations for ionic propa- 
gators, self-energies and vertex functions includes many 
such integral equations. Great care is needed to identify 
all different processes together with the proper assign- 
ment of quantum numbers, energies, matrix elements and 
signs. This coupled system is solved self-consistently for 
the unknown self-energies and ionic propagators. There 
are between 6 and 158 different ionic propagators, de- 
pending on symmetries, and about 10-100 different dia- 
grammatic contributions to each self-energy. Therefore, 
the number of necessary convolutions in each iteration 
step is easily in the thousands. As in the usual proce- 
dure, an additional set of coupled integral equations has 
to be solved for the so-called defect propagators. ^^'^^ 

While the computational effort for the two-impurity 
schemes presented here is increased considerably com- 
pared to their single-impurity counterparts, the numeri- 
cal difficulties are of the same nature: The ionic propaga- 
tors exhibit very narrow features near the ionic thresholds 
which can be handled by the usage of optimally adapted 
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FIG. 2. Contributions to the self-energy Eoo;Oo(2)- The squares represent ionic propagators, Pmjimj2i'nii»"i2 (•^)' where the 
earlier [later] state characterized by the quantum numbers Af — {mii,mi2) [M = (myi, 771/2)] is indicated by the attached lower 
[upper] two lines. For clarity, a diagonal line inside the box is added in case of a non-diagonal propagator. Summations over 
all internal quantum numbers is understood in all graphs presented in this work. 
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FIG. 3. Contributions to the self-energy EctO;Oo-(2). The double-line visible in some contributions indicates a doubly occupied 
ionic shell on one site. 



energy meshes. The strongest challenge is posed by the 
numerical solution of the system of equations for the de- 
fect propagators. The lack of an absolute scale, i.e. the 
lack of the knowledge of the partition function, leads to 
suboptimal convergence properties. ^^'^^ However, proper 
and rather accurate solutions can indeed be obtained, 
and results will be presented in the following sections. 
The interested reader can consult Refs. 25 and 27 for 
more details. 

When extending these approximations to systems with 
more local orbitals or impurities one is confronted with 
the fundamental problem of an exponentially increasing 
number of local many-body states. Additionally, the 
inclusion of vertex corrections as in the two-impurity 
ENCA further increases the contributions to each ionic 
self-energy. However, in the present approach symme- 
tries can easily be implemented by identifying equivalent 
propagators; thus, the computational cost is dramatically 



reduced. In this respect it is also very effective and eas- 
ily implemented to restrict the allowed local occupation 
and to exclude many-body states with very high excita- 
tion energies, as it was done for the two-orbital Anderson 
model (Ref. 22). 



IV. RESULTS FOR UNCOUPLED IMPURITIES 

In this section we first concentrate on the basic ques- 
tion of how the novel two-impurity solver performs in 
some simple model situations. We compare the two- 
impurity SNCA and ENCA with the corresponding 
SNCA- and ENCA-solutions of the single-impurity An- 
derson problem in order to check on the presence and 
relative weight of artificial correlations introduced by the 
two-impurity treatment. 

We assume spin-degeneracy in the following and mea- 
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FIG. 4. Contributions to the vertex Aoo;ctO within the two-impurity ENCA approximation which includes the first oder of vertex 
corrections. 



sure energies with reference to the chemical potential 
fi = uj = 0. We always take the noninteracting con- 
duction electrons to have the tight-binding dispersion of 
a three-dimensional simple-cubic lattice 



-'2t[cos{kx) + cos(fcy) 4- cos{kz)] 



(22) 



(with the lattice spacing set to one, a = 1). As the unit 
of energy we take either the half-bandwidth of the non- 
interacting conduction band, D = 6i, or the Anderson 
width A, depending on what seems more appropriate. 
We also restrict ourselves to identical impurities, that is 

U, and = pj and mostly 



focus on the two-electron sector. 

We define even (+) and odd ( 
states^" 

1 



parity one-particle 



V2 



(23) 



and also introduce even and odd hybridization functions, 
^+(^) = I E m^^os" {\kd)i,5{u^ ~ el) (24) 

k 



.1 



N 



Y,m^s\n\-kd)^5{u^el). (25) 



Here, d^Ri — R^ is the distance vector between the two 
impurities. For very large distances, \d\ — >■ oo, IS.±{ijj) 
both approach the single-impurity hybridization func- 
tion. 



(26) 



For vanishing distance, \d\ 0, A+(w) = 2A'=^^^(i:j) 
and A_(a;) = 0. 

These even and odd hybridization functions arc shown 
in Fig. 5 for the three-dimensional tight-binding band- 
structure of a simple-cubic lattice and the distance vec- 
tor along a principal direction, e.g. d = de^. We also as- 
sumed a constant hybridization matrix element, Vk = V. 

The first limiting case we like to investigate is obtained 
if all couplings between the impurities of the cluster are 
suppressed, i.e. A12 — 0, and likewise all matrix elements 
in the Hamiltonian which directly couple the two impu- 
rities, i.e. U{2 = and = 0. 

Although the resonant level limit [/ = is trivially 
solved exactly, the reconstruction of this solution via di- 
rect perturbation theory is highly nontrivial.^^'^^ In fact, 
this case turns out to be particularly unfavorable for di- 
rect perturbation theory. 

The spectral densities of even and odd states in the 
case of [/ = can be easily derived from the exact one- 
particle Green functions 



G±iz) 



1 



lmG±{u + iS) 



1 



where generally an additional coupling function 



9{z) 



1 \Vk?suv{kd) 



N 



(27) 
(28) 

(29) 



occurs. However, for the present case of a momentum 
independent hybridization matrix element, Vk = V , and 
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inversion symmetric lattices where 
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FIG. 5. (Color online) Hybridization functions normalized to 
the Anderson width A for the even (a) and odd (b) channel for 
various impurity-impurity distances d = |d| and momentum 
independent hybridization matrix element, \Vk}^ = \V\^ . For 
comparison the single impurity hybridization function is also 
shown, which is nothing but the rescaled density of states of 
the noninteracting band electrons. 

The exact curves for U ~ 0, = —OAD and distance 
d ~ \d\ ^ 1 are compared with the results of the two- 
impurity SNCA and two-impurity ENCA calculations in 
Fig. 6. Although the overall features of all curves show 
rough agreement, the two approximations perform dif- 
ferently in detail. Both capture the van-Hove singularity 
on the right flank of the resonance at ui/D = — i, but 
the two-impurity ENCA comes closer to the exact solu- 
tion than two- impurity SNCA. The position and height of 
the maximum is also captured better by the two- impurity 
ENCA. Additionally, the two-impurity SNCA produces 
some unphysical weight in the left flank of the peak at 
Lj/D w 0.75. 

As mentioned in the previous section such deficiencies 
had to be expected. They turn out, however, to be less 
pronounced as one might have feared, so that even in this 
most critical case of uncorrelated impurities a reason- 
able result is produced by the two-impurity ENCA. As 
in the single impurity case, the performance should even 
increase when correlations become important for U > 0. 
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FIG. 6. (Color online) Exact spectral functions of the res- 
onant level with — —OAD [see Eq. (27)] in the even (a) 
and odd (b) channel for distance d — 1 compared to the two- 
impurity (TI-) SNCA and ENCA resuhs. 



We now turn to the Kondo regime with U > — 2e-' > A. 
The hybridization V tends to produce separate Kondo ef- 
fects at each impurity controlled by the diagonal elements 
of the hybridization functions Ajj^ui), but additionally, 
the non-diagonal element of the hybridization, A;j(cj) 
(with I j), will produce a RKKY- interaction between 
the local magnetic moments on the impurities. In order 
to benchmark the approximations we neglected the latter 
for the moment and set Aij{uj) = for I ^ j- 

Calculations for the Kondo regime arc presented in 
Figs. 7(a) and 7(b). Care has to be taken to facilitate 
a comparison with the results of single-ion SNCA- and 
ENCA-solvers, which are shown in these figures, too. 
Since our two-impurity SNCA and ENCA calculations 
treat the two-impurity systems, though not coupled, as a 
whole, the groundstate energy and threshold energies for 
local excitations should have twice their single-ion value, 
and the single-particle spectral functions are normalized 
to two. Hence the energies and curves are correspond- 
ingly scaled with a factor of two to be comparable. 

Figure 7(a) shows the spectrum of the ionic propagator 
in the energy region around the ionic threshold for singly 
occupied impurities, i.e. ppi{uj) = ~—lm[Pa,„{uj + iS)] 
for the SIAM and ppi{oj) = -ilm[pL',^^' (w/2 + id)] 
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FIG. 7. (Color online) (a) Threshold behavior of the differ- 
ent SIAM and the decoupled TIAM approximations deduced 
from the ionic propagator for one electron per site and for 



■f , [/ = D, T = 0.003D, and A = O.ID. For the 
SIAM ppi{oj) = - j;lm[Pcr.^{oj + iS)] and for the TIAM as 
pPi{u}) — —■^lm[Pcrcr',iTa-'{^/'2 + iS)]. (b) FuU one-particle 
spectral functions of the different approximations for the same 
parameters as in (a). The inset shows a close-up of the region 
around the Fermi level. 



onances of the TIAM are less pronounced. All curves 
are calculated with the same absolute temperature T = 
0.0003J3 which implies different relative temperatures to 
the corresponding Kondo scale of each approximation. 
The TIAM solutions produce a slightly smaller Kondo 
scale and consequently the peaks are not as high, which 
confirms the conclusions drawn from the ionic propaga- 
tors. 

The good qualitative and mostly also quantitative 
agreement of the calculations support our confidence in 
the validity of the results presented for coupled impuri- 
ties in the next section. 



V. RESULTS FOR COUPLED IMPURITIES 

The physics of a fully coupled cluster of two impuri- 
ties can become amazingly rich when all possible sources 
of interactions are taken into account. Wc therefore re- 
strict ourselves to a few important cases, which we first 
will consider separately in order to clarify basic physical 
effects. 

In general, the hybridization part of the Hamiltonian 
will not only induce the tendency towards a local Kondo- 
scrccning of moments considered in the previous sec- 
tion, but will also induce an additional effective RKKY- 
cxchange interaction with the asymptotic form in leading 
order perturbation theory^* 

^s'm{2kFd) + 2kFdcos{2kFd) 

JRKKYOC-V/ 4 , (jUj 

(z/cpctj 

where hp is the Fermi wave vector along the distance 
d 3> 1 between the impurities (lattice constant a = 1). 
The RKKY-interaction according to Eq. (30) oscillates 
with distance and can induce a ferromagnetic or antifer- 
romagnctic coupling. 



for the two-impurity model. As is well known^^ the 
ENCA performs considerably better than the SNCA 
since, for instance, it produces a more accurate, i.e. 
larger, Kondo temperature Tk- Tk can be estimated 
as the difference between the threshold, i.e. the position 
of the maximum in ppi, and the single-particle energy, 
Tk ^ -^threshold ~ c^- The same apparently holds when 
comparing of the two-impurity SNCA and two-impurity 
ENCA approximation. The two-impurity ENCA even 
seems to exceed the single-impurity SNCA and comes 
near to the quality of the single-impurity ENCA, at least 
for the typical parameter values shown. 

A similar conclusion can be drawn from the one- 
particle spectral functions of Eq. (13) shown in Fig. 7(b). 
They exhibit the typical three-peak structure with a 
Kondo-resonance near the Fermi-level and large spectral 
weight at the so-called Hubbard peaks which originate 
from the ionic one-particle energies e-^ and e-^ + U . While 
the high-energy features of the SIAM and TIAM solu- 
tions are in good agreement, the low-energy Kondo res- 



A. Isolated two-impurity cluster 

It is instructive to consider, what physics already is 
included without hybridization V to the bandstates. For 
this purpose one can inspect the exact solution of the 
eigenvalue problem for the isolated cluster (V^ = 0) in 
the presence of interactions which amounts to setting all 
hybridization functions to zero, Aji{uj) = 0. We include a 
direct magnetic exchange term and also allow for a direct 
single-electron transfer, 

= E *(/l/2. + h.c.) ~ Jsisi. (31) 

Since the exchange interaction has no effect on the one- 
particle sector of the cluster, the one-particle eigenstates 
correspond to the (anti-)symmctrized operators given in 
Eq. (23) with eigenenergies 

=ei = e^ ±t. (32) 
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The two-particle sector contains six states, three of 
which belong to a degenerate triplet with one electron 
on each of the sites. For these states spin-conserving 
hopping is ineffective due to the Pauli-principle, which 
leads to: 



1,771 = cr) |cr)i|cr)2, 



(33) 



|5=l,,n = 0) = 4sflt>lU>2 + U)l|t>2 



with eigenenergies 



The three remaining singlets in the two-particle sub- 
space are 



IS" = 0,777 = 0)o = 



^(it,;>i|o)2 



L|t,;>2 (34) 



E 



(2) _ 



So 



2e^ + U, 



and 



|5 = 0,777 = 0)± = — 



1 

N± 



t)ii;)2-i;)iit): 



+ a± it,;)i|o)2 + |o)iit,;); 



E 



(2) 

S± 



3J 



2e + -r - «± 



(35) 



where 




N± 



In order to understand the result it is useful to realize 
that for the polar state \S = 0,r7i = 0)o direct transfer 
is irrelevant due to orbital antisymmetry. The effect of 
direct hopping on the remaining two singlets, as well as 
the influence of direct exchange J are most clearly seen 
in the limiting case where U exceeds t and J, i.e. U ^ 
\t\, I J|, where one finds 



i(2) 

s+ 

i(2) 



4 U 



2ef 



(36) 
(37) 



In this limit the well known competition between the 
triplet of Eq. (33) and the singlet of Eq. (37) for becoming 
the ground state is realized. The physics is determined 
by the exchange splitting of these levels. 



AE ^ E 



(2) 



E, 



(2) 



4t2 

J- — = J + Jt. (38) 



For larger values of the direct transfer t the level- 
splitting becomes increasingly important and renders a 
discussion based exclusively on an effective exchange in- 
teraction impossible. 

States with three or four electrons in the cluster have 
higher energies through additional contributions of U. 
They can be deduced from the above results by particle- 
hole transformation. 



B. Direct exchange coupling 

For the following numerical study we will set t = but 
consider a direct coupling J. We also include a finite di- 
agonal hybridization but still ignore the non-diagonal hy- 
bridization functions producing the RKKY interaction, 
i.e. Ai2(a;) = 0. Therefore, we choose a direct two- 
impurity interaction Hamiltonian [see Eq. (2)] 



-f-f 



(39) 



Then, the typical competition between the single-ion 
Kondo effects and the fixed nonlocal exchange J results. 
The Kondo coupling between each impurity spin and the 
conduction electrons is governed by the antiferromag- 
netic exchange 



7 2[/A 

"-'Kondo f/f,rr\'^^ 

eJ (eJ + U) 
leading to a characteristic energy scale 



K 



kondo CXp ( — Tt/ JKondo) 



(40) 



(41) 



with a w 27rmin{J7, D}. 

Onc-particle spectral functions calculated with the 
two-impurity ENCA-solver are shown in Fig. 8 for vari- 
ous values of a ferromagnetic direct coupling ( J > 0) and 
temperatures. The overall form of the spectral function 
and the high-energy features are essentially unchanged 
under the inclusion of J. However, pronounced changes 
occur in the low-energy region around the Fermi energy 
cj = 0, which at J = and small T < Tk exhibits the 
well-known many-body resonance of the single-ion Kondo 
effect. 

For very large direct exchange coupling J > the two 
impurity spins are aligned parallel and nearly act as one 
rigid spin with S = I. This corresponds to the triplet 
state in the isolated cluster which experiences a modified 
two-stage nonlocal Kondo screening with a reduced char- 
acteristic energy scale y^-^'™?' < Tk.^^ This leads to the 
formation of a narrower and less saturated Kondo reso- 
nance at w = 0; clearly visible in Fig. 8(b) and is similar 
to what happens in multi-orbital SIAM.^^ 

At excitation energies of the order of the exchange 
coupling, u! « ±J, one probes a domain dominated by 
correlations typical for the single-ion Kondo effect. The 
weaker nonlocal correlations of the low-energy regime are 
then broken. The two impurities appear as essentially 
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FIG. 8. (Color online) Spectral function for e-' = -f, (7 = D, 
A = 0.1-D with direct ferromagnetic coupling J > 0: (a) 
Fixed temperature T = 0.004D = 0.04A and various values 
of the coupling, (b) Magnification of the region around the 
Fermi-energy, (c) The spectral function near the Fermi level 
for fixed coupling J = 0.017-D = 0.17A and various tempera- 
tures. Notice, that while in panel (a) the energy is measured 
in units of the half-bandwidth D, panels (b) and (c) use the 
Anderson width A. 



uncoupled and behave like isolated single impurities. In 
the spectral functions this is indicated by the humps vis- 
ible for small J which develop into side-maxima at larger 
J vifhich are located the energies w « ± J. For energies 
larger than > J the J = spectral function is ap- 
proached. This is in accord with the finding, that the 



characteristic high-energy scale of this model is essen- 
tially given by the single-impurity Kondo scale. 

This is also reflected in the temperature dependency 
as shown in Fig. 8(c). The many-body resonance at first 
forms like in a SI AM, and only for low temperatures T < 
J the narrowing near the Fermi level sets in and causes 
the above mentioned humps or maxima. 
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FIG. 9. (Color online) (a) Spectral functions for various an- 
tiferromagnetic direct exchange couplings J < and fixed 
temperature T = 0.04A. Panel (b) shows a magnification 
of the low energy region, (c) The spectral function around 
the Fermi level for fixed coupling J = — 0.17A and various 
temperatures. Notice, that while in panel (a) the energy is 
measured in units of the half-bandwidth D, panels (b) and 
(c) use the Anderson width A. Other parameters are as in 
Fig. 8. 
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We now turn to the case of an antiferromagnetic direct 
coupling (J < 0) where the nature of the physical sce- 
nario for low lying states changes qualitatively. Now a 
real competition takes place between local singlet forma- 
tion via the Kondo-effect and nonlocal singlet binding in 
the molecular two-electron state. The single-ion Kondo 
effect present at small J is suppressed with increasing 
negative values of J . This is clearly borne out by the 
spectral functions of Fig. 9, which exhibit a very rapid 
depletion of the spectral function near the Fermi energy 
with |J|. 

The formation of a gap around w = then implies that 
the low lying molecular singlet does effectively not inter- 
act with the band states. However, side peaks at the 
edges of the gap can be observed at energies w = ±J. 
These are indicative for the fact that the excited molec- 
ular triplet experiences a kind of virtual Kondo effect.'^" 

This interpretation is supported by the temperature 
dependence of the curves depicted in Fig. 9(c) for a fixed 
value of J = —0.1 7 A. In the same way as the pseudogap 
around a; = forms for temperatures lower than | J|, the 
side peaks emerge and increase in height. 

C. Direct single-particle hopping 

We now neglect the direct magnetic exchange, but in- 
troduce a direct single-particle hopping, i.e. 

Wf ^tY^Ulha + h.c). (42) 

We still ignore the non-diagonal hybridization 

Ai2(w)=0. (43) 

Based on the discussion of the isolated two-impurity 
cluster in Sect. VA we can expect a twofold source for 
modifications of a pure single-ion Kondo effect: (1) t gives 
rise to an effective antiferromagnetic exchange interac- 
tion Jf = — and (2) t splits the ionic one-particle 

levels — >■ e£ and produces a tendency towards even 
and odd molecular onc-particlc states. Therefore, the in- 
teresting question arises, how the scenario developed for 
an antiferromagnetic direct exchange interaction will be 
altered by the impending effect of even-odd splitting. 

Figure 10 shows the spectral function for various val- 
ues of t. One recognizes the expected opening of the 
pseudogap around the Fermi level w = with increasing 
absolute value of Jt [see panel (b)]. 

The hopping induces changes of the spectral function 
at all energies [see 10(a)], which is in contrast to the case 
of a direct coupling where only the low-energy region 
was affected. Particularly, and in accord with our expec- 
tations, the ionic Hubbard resonances near oj = and 
w = + U are split by the hopping. The separation 
of the peak maxima is roughly given by 2t which is in 
accord with the simple reasoning from the isolated two- 
impurity model, see Eq. (32). This reflects a tendency. 
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FIG. 10. (Color online) (a) Spectral functions at a fixed tem- 
perature T = 0.04A for various values of the direct single- 
particle hopping t. To reveal its physical influence more 
clearly the effective antiferromagnetic exchange Jt = —^J- is 
given instead of the bare hopping t. For example, Jt = — 0.4A 
corresponds to t = O.ID and Jt = —0.7 A corresponds to 
t = 0.13Z). Panel (b) shows a magnification of the low energy 
region (energy is measured in units of the Anderson width 
A). Other parameters are as in Fig. 8. 

which has been outlined before''^, namely that features 
of the unperturbed structures of the one-particle states 
leave their traces in the quasiparticle bandstructure. 

In the low-energy region, hopping produces a gap simi- 
lar to the case of a direct exchange coupling [see Figs. 9(b) 
and 10(b)]. However, whereas the gap is not as pro- 
nounced as in Fig. 9(b), the side peaks at the edge of the 
gap are considerably higher. 

The question arises whether the splitting is due to a 
suppression of the Kondo effect as in the case of a direct 
antiferromagnetic J. It can be addressed by decompos- 
ing the spectral function into even and odd parity parts 
as would be produced by the operators /J_^ of Eq. (23). 
This is shown in Fig. 11 for t = O.ID, i.e. Jt = -0.4A. 
All the peaks in the spectral function can clearly be at- 
tributed to (mostly) one parity channel. In particular, 
the quasiparticle Kondo resonance is partitioned into an 
odd contribution for negative energies and an even con- 
tribution for positive energies. 

This suggests that the splitting is not exclusively due to 
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FIG. 11. (Color online) Spectral function of Fig. 10(a) for 
Jt/A = —0.4 decomposed onto even and odd contributions. 



a suppressed Kondo effect, but is mainly a single-particle 
even-odd splitting of the quasiparticlc excitations. The 
side peaks then indicate the local preformation of quasi- 
particlc bands, with the odd states associated with the 
bonding region at A; = and the even states with the 
anti-bonding region at /c = ±7r. 
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FIG. 12. (Color online) (a) Low-energy gap as function of 
\Jt\. The parameters for the dashed fit are tc — 0.05 A and 
a — 2.16A (b) Occupation numbers of the singlet and triplet 
states in the two-particle sector as function of the (negative) 
effective antiferromagnetic exchange Jt- All other parameters 
are as in Fig. 10. 

We show the width of the central low-energy pseudo- 
gap as function of the effective exchange \Jt\ = ^ in 
Fig. 12(a). The gap first opens at a finite value tc for 
the hopping. For smaller values of Jt the splitting be- 
haves like a square root indicating an opening of the gap 
which is linear in t. At larger t the gap is dominated by 
the impact of the quadratic i-dcpcndence of the exchange 
coupling and thus goes linear in \Jt\. Indeed a fit with 

a function a{t — tc) + 4 ^ jj"' reproduces the gap quite 
well. 

The molecular triplet and singlet states are separated 
in energy roughly by \Jt\ [see Eq. (38)]. At large t the 
singlet becomes lower in energy that the triplet. This 
can be directly observed in Fig. 12(b), where the oc- 



cupation numbers of these molecular states are shown 
as function of |J(|. For small \Jt\ the occupation num- 
bers are roughly equal and mixed at a nearly constant 
ratio. With increasing coupling, the singlet occupation 
rises steeply, accompanied by a corresponding fast de- 
pletion of the triplet states. The large- |Jt| regime thus 
reproduces the behavior already known from a direct an- 
tiferromagnetic exchange. 

D. Non-diagonal hybridization and RKKY 
exchange 

It is left now to clarify the role played by the non- 
diagonal hybridization. We set the direct parameters t 
and J to zero and include the full hybridization function 
of Eq. (20), which allows for a coherent hybridization of 
the cluster states with the conduction band. The non- 
diagonal component for distance d = i?2 — i?i reads 

Ai2(w) =^ cos{kd)TTS{u - el). (44) 

k 

The process of non-diagonal hybridization is qualitatively 
different from a direct hopping t as considered before. As 
a consequence of the intermediate propagation through 
the band, a nontrivial dependence on the excitation en- 
ergy Lu results. The hybridization function exhibits an in- 
creasingly oscillatory behavior with growing distance be- 
tween the impurities while the absolute height decreases 
as shown in Fig. 13. 
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FIG. 13. (Color online) Non-diagonal hybridization function 
of Eq. (44) for various distances as indicated. The inset 
sketches the three-dimensional free density of states of the 
conduction electrons. 

The following effects are expected as consequences of 
diagonal and this non-diagonal hybridization: (1) the 
RKKY-exchange interaction between the two impurity 
spins is induced, which favors ferro- or antiferromag- 
netic alignment depending on the distance. (2) Kondo 
screening-clouds (see, for example, Ref. 32) for singly oc- 
cupied impurities are induced in the band. In case the im- 
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purities are spatially not well separated and the individ- 
ual Kondo-clouds penetrate each other, this then involves 
considerable nonlocal coherence. (3) A tendency toward 
the stabilization of even and odd parts of the molecular 
single-electron orbitals might appear, as the result of an 
induced (indirect) transfer via the band states. 

Thus, we do not expect a simple realization of the 
Doniach-scenario of competing local and antiferromag- 
netic correlations, which would either lead to a Fermi 
liquid with quenched spins (and some short-ranged anti- 
fcrromagnetic correlations) or to antiferromagnetic order 
between the impurities. 
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FIG. 14. (Color online) Impurity spectral functions with di- 
agonal and non-diagonal hybridization functions for two dis- 
tances (a) |d| = 1 and (b) |d| = 2 for various Anderson widths 
as indicated. The inset shows a close-up of the low energy re- 
gion. Other parameters are: = U = D, T = 0.004D 
and J = t = 0. 

In Figure 14 the one-particle impurity spectral func- 
tions for two different impurity distances along the same 
principal direction of a three-dimensional simple-cubic 
lattice are shown for various values of the hybridization 
strength A. These cases can loosely be associated with 
situations already studied above. Our choice of parame- 
ters places the nearest-neighbor d = 1 case into a regime 
of antiferromagnetic exchange, Jrkky < 0, whereas at 
d = 2 the impurities should experience a ferromagnetic 
^RKKY > 0. 

The spectra greatly change with increasing hybridiza- 



tion. The Hubbard peaks are broadened and are moved 
to larger energies, as it is expected. These high-energy 
features are mainly determined by the diagonal part of 
the hybridization and are insensitive to the distance be- 
tween the impurities, thus the figures of panel (a) and 
(b) are very similar at large |a;|. Increasing A first im- 
plies a larger absolute value of Jrkky < for d = 1. 
Consequently, a pscudogap forms at the Fermi level for 
A > 0.151? as visible in the inset of Fig. 14(a). The rea- 
son for the gap being not visible at smaller A is the finite 
temperature of T = 0.004Z? used in the calculation. 

For distance d = 2 the larger ferromagnetic Jrkky 
should lead to a narrowing of the many-body reso- 
nance with side-peaks, but this is not clearly observed 
in Fig. 14(b). On the other hand, the local single- 
impurity Kondo scale is also enhanced with increasing 
A [see Eq. 41], leading to an increase in the overall spec- 
tral weight of the Kondo resonance at the Fermi level. 
(This is also observed in the case of the antiferromag- 
netic d = 1 situation.) The signatures of indirect fer- 
romagnetic exchange are by far not as clear as in the 
case studied in a previous section with direct exchange. 
For one, they are strongly competing with the increas- 
ing single-ion Kondo effect and for another, the dynamic 
non-diagonal exchange does not act in the same simple 
way as an static exchange constant. Therefore, the Do- 
niach scenario seems to be too oversimplifying. 
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FIG. 15. (Color online) Impurity spectral functions decom- 
posed into even [(a) and (c)] and odd [(b) and (d)] contribu- 
tions for distance |d| = 1 [(a) and (b)] and |d| = 2 [(c) and (d)] 
for various Anderson widths. The inset shows a close-up of 
the low energy region. All other parameters are as in Fig. 14. 

The influence of molecular correlations favoring even 
and odd one-particle states via an induced effective trans- 
fer can again be derived from a decomposition of the spec- 
tra in the even and odd parity components. These are 
shown in Fig. 15 for the same parameters as in Fig. 14. 
It is obvious, that the indirect hopping is much less effec- 
tive in producing even-odd splitting when compared to 
the direct hopping of the previous section. The even and 
odd spectral functions are very similar to each other and 
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only for the small distance d = 1 moderate asymmetries 
occur. The splitting is noticeable, though not nearly as 
strong as in Fig. 11. But interestingly, the splitting of the 
low-energy quasiparticle peak at the Fermi level is much 
more pronounced than the imbalance of the high-energy 
Hubbard peaks. This supports the idea that nonlocal 
correlations cause the coherence in the low-energy region 
where the Kondo effect develops. They also lead to pseu- 
dogap formation and in consequence to the notion of a 
impending quasiparticle bandstructure, a viewpoint al- 
ready adopted above. 

For d = 2 such an effect is barely visible and an asym- 
metry can not be observed in the spectral functions. The 
odd channel has a slightly larger spectral weight in the 
many-body resonance around the Fermi level as com- 
pared to the even spectra. 




FIG. 16. Height of the two-impurity spectral function at 
uj — relative to its value without the second impurity [see 
Eq. (45)]. (a) A = O.UD (b) A = 0.16D. One impurity is 
located at the origin, while the other at the distance vector 
d = djje^ + dyEy. The insets show the sign of this quantity, 
where black denotes Ap < and white Ap > 0. All other 
parameters are as in Fig. 14. 

The investigation of the effective transfer- and ex- 
change can easily be extended to other relative positions 
of the two impurities in a cubic lattice. As one might 
expect, an oscillation of physical quantities which can be 



associated with ferromagnetic or antiferromagnetic ten- 
dencies is found with varying distance. One example of 
such an oscillating quantity is the height of the many- 
body resonance right at the Fermi level in comparison to 
its value without the second impurity, i.e. 

dp = p{LJ = 0)^P\d\=ooiL^^O). (45) 

We present this quantity in Fig. 16 for two different hy- 
bridization strengths as a function of the two-impurity 
distance d = dxG^ + dyCy in a two-dimensional plane, 
where one impurity is always located in the origin. Red 
and yellow colors imply a larger height of the two- 
impurity spectral function than in a SIAM, while black 
and blue colors indicate a suppression of p(uj = 0). 

Regions with enhancement and suppression of the 
many-body resonance can be clearly identified. While a 
suppression can easily be explained by the antiferromag- 
netic exchange coupling, the reason for an enhancement is 
not as clear. We attribute it to the effect of the dynamic 
non-diagonal hybridization, producing a constructive in- 
terference and an effectively larger total hybridization. 
However, the spatial pattern is not simply determined 
by the Manhattan distance, i.e. the number of elemen- 
tary hoppings between the two impurities. It appears, 
that more complicated interferences of single-particle and 
interaction effects obviously play a role. For very large 
distances these short-distance effects should be irrelevant 
and an asymptotic form like Eq. (30) is expected. 

VI. CONCLUSION 

In this work we have presented a novel two-impurity 
solver which is based on direct perturbation theory with 
respect to the hybridization. It is capable of treating any 
kind of direct interaction of the two-impurities as well as 
the dynamical indirect coupling via the conduction band. 

Including separately different couplings of the two im- 
purities we carefully investigated the role played by each 
and the resulting physical mechanisms. In situations, 
where the impurities are not coupled via the conduction 
electrons but only directly, we found that a ferromag- 
netic exchange leads to a narrowing of the many-body 
resonance which is characteristic of a higher-spin Kondo 
effect. 

For antiferromagnetic coupling, the inter-impurity sin- 
glet formation competes with two separate local Kondo 
effects. This leads to the suppression of the Kondo effect 
due to singlet-triplet splitting at small enough tempera- 
ture which is clearly signalled by a pseudo-gap opening 
in the many-body resonance at the Fermi level. 

Interestingly, such a pseudo-gap is also induced with- 
out direct exchange coupling but instead including a di- 
rect single-electron hopping between the impurities. This 
produces a splitting of bonding and anti-bonding two- 
impurity orbitals, i.e. of even and odd parity components 
of the spectral function. Additionally, an effective anti- 
ferromagnetic exchange is generated by a finite hopping. 
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This also supports the formation of a gap in the many- 
body resonance which, however, is accompanied by sharp 
side-peaks of the even and odd excitations. 

In the generic scenario without direct interactions be- 
tween the impurities but with couphngs generated indi- 
rectly via the conduction electrons the situation is not as 
clear. Signatures of the induced RKKY interaction are 
observable in form of oscillatory behavior with varying 
distance between the impurities. It seems that the in- 
duced magnetic exchange interactions are the dominant 
source of competition with isolated Kondo effects. Molec- 
ular correlations and splitting of even and odd parity 
channels induced by an effective single-particle hopping 
between the impurities are only observable for smallest 
distances. Even though the total size of these effects is 
rather small, the low-energy many-body physics seems to 
be very sensitive to this kind of perturbation. 

Although the well known Doniach-scenario captures 
essential aspects of the physics involved, we are led to 
the conclusion that the coupled two-site cluster exhibits 



a much richer physical behavior, in particular at small 
impurity distances. For this reason the dynamical non- 
diagonal hybridization function needs to taken into ac- 
count. 

The ability of the solver to work with arbitrary two- 
impurity distances and arbitrary hybridization functions 
allows for utilization in a nonlocal extension of dy- 
namical mean- field theory. '^"^ We propose such a scheme 
where two-site correlations are included in a separate 
publication'^* and the results with the two-impurity 
ENCA as impurity solver are very promising. 
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